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Abstract 

. The study of social networks and especially of the stochastic dynamics of the diseases spread in 

human population has recently attracted considerable attention in statistical physics. In this work 
we present a new statistical method of analyzing the spread of epidemic processes of grippe and 
^ ■ acute respiratory track infections (ARTI) by means of the theory of discrete non-Markov stochastic 

processes. We use the results of our last theory (Phys. Rev. E 65, 046107 (2002)) to study statistical 
memory effects, long - range correlation and discreteness in real data series, describing the epidemic 
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CD ■ dynamics of human ARTI infections and grippe. We have carried out the comparative analysis of 

'nT ■ 

I the data of the two infections (grippe and ARTI) in one of the industrial districts of Kazan, one 

of the largest cities of Russia. The experimental data are analyzed by the power spectra of the 
^ ' initial time correlation function and the memory functions of junior orders, the phase portraits of 

the four first dynamic variables, the three first points of the statistical non-Markov parameter and 
the locally averaged kinetic and relaxation parameters. The received results give an opportunity 
to provide strict quantitative description of the regular and stochastic components in epidemic 
dynamics of social networks taking into account their time discreteness and effects of statistical 
memory. They also allow to reveal the degree of randomness and predictability of the real epidemic 
process in the specific social network. 
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The functioning of network-organized statistical systems essentially depends 
on the nature of interaction between their elements. It is especially due to the ef- 
fect of disease-causing contacts and topology of networks. This reason as well as 
the variety of the displayed nonlinear behavior have made this problem the sub- 
ject of several studies by fundamental methods of statistical physics. To analyse 
the epidemic and disease dynamics complexity, it is necessary to understand the 
basic principles and notions of its distribution in long-time memory social media. 
Here we consider the problem from a theoretical and practical viewpoint and 
present the quantitative evidence confirming the existence of stochastic long- 
range memory and robust chaos in a real time series of respiratory infections 
of human upper respiratory track. We also discuss the implications of discrete 
non-Markov stochastic processes in real social complex systems from the point 
of view of the recent theory (Phys. Rev. E 65, 046107 (2002)). 



I. INTRODUCTION 



In the last few years it has been well recognized that the study of epidemic and disease 
dynamics in social networks is a relevant theoretical issue concerning the spreading of viruses 
and diseases of various nature It is connected with the phase transitions in 

models of agent spreading in nonsteady systems. In particular, it concerns the development 
of the models describing an epidemic spread, forest fires, growth of populations, activity of 
catalyzers, formation of stars and galaxies. The given problem might be of interest for all 
those interested in physical problems, related to the critical phenomenon in the universality 
class of Reggeon field theory, contact processes, ordinary and directed percolation and scale- 
free properties in many real systems such as Internet, World Wide Web, food webs, protein 
and neural networks. 

This study is of significant interest for modern biomedicine, ecology and economy 



22j | . The epidemic means that a large number of people in a certain country simultaneously 



develop a disease at a certain period of time. In case of "pandemic" people in several 
countries develop the same disease at the same period of time. The first cases of respiratory 
infections with fatal outcome in the USA and Europe were described at the beginning of the 
20th century. About 20 million people died then of " Spanish fever" within the period of 10 
months. The pandemic of "Asian" and "Hongkong" grippe with smaller quantity of fatal 
outcomes [^-0 were observed in the middle of the 20th century. The period of pandemic 
processes is about 30-40 years. As a rule, epidemics break out in the autumn or in the winter 
in Northern Hemisphere and in the spring or in the summer in Southern Hemisphere. The 
duration of epidemics is about 1-3 months. Seasonality is one of the typical demonstrations 
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of the phase development of epidemic processes. Respiratory infections affect people of any 
age. However children aged 1-14 develop the disease four times more often than adults. Each 
epidemic process of respiratory infections is accompanied by the increasing death rate. The 
parameters of grippe death rate in the world is 0, 01 — 0, 2%. The greatest number of deaths 
caused by grippe is connected with serious complications after the infection. Children (up 
to 2 years old) and elderly people (over 65 years old) are more liable to this infection and 
die more often. 

Various authors have discovered that the number of people developing the disease de- 
pends to a great extent on ecological and economic conditions. The most common diseases 
the population of the earth suffers from are infections of respiratory organs among which 
grippe and acute infections of the upper respiratory track are main components. The low 
standard of living as well as the pollution of the environment have had their negative in- 
fluence on people's health. This is the main reason for the worsening epidemic situation in 
the world. Respiratory infections cause significant economic damage. Expenses caused by 
epidemics affect successful development of economy due to the expenses needed for financ- 
ing scientific actions to prevent epidemics. These expenses are commensurable only with 
the sums earmarked to prevent heart diseases In many countries of the world great 

money is spent on the development of vaccines and serums, as well as prevention of epidemic 
outbreaks. However the viruses of grippe and some other respiratory infections have great 
variability. There is a great variety of grippe virus cultures. Therefore the coming grippe 
variant is difficult to predict. The identification of the virus presents a great problem, conse- 
quently the prevention of epidemic outbreaks is complicated. Thus new preventive measures 
and methods of predicting possible epidemic outbreaks 17|- 2^ are of great interest. The 
unexpected outbreak of grippe in March 2003 is a good example of the case. 

Now various methods of description and prediction of epidemic and pandemic processes 
are developed in medicine [l3-22|. Hypotheses about risk factors have been formulated 
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as a result of these studies. The standard statistics indices 
epidemic studies. The basic data yield the number of diseases, their intensive parameters 
and average sizes. The extensive parameters, the cumulative data, the relative number of 
presentation, the ratio parameters and standardized parameters are used depending on the 
specific features of the epidemic process. However existing methods are not efficient enough. 
A model for the spread of an infection is analyzed [2^ for different population structures. 
For more ordered systems, there exists fluctuating endemic state of low infection. At a finite 
value of the disorder of the network, a transition to self-sustained oscillations in the size of 
the infected subpopulation has been recorded. A spatial model related to bond percolation 
for the spread of a disease that includes variation in the susceptibility to infection has also 
been considered in Ref. |24|. 
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Methods of networks are all important for fundamental solution of the problems related to 
the spread of epidemics of grippe and ARTI. Traditionally these systems have been modelled 
as random graphs, however topology and evolution of real networks are governed by robust 
organization principles. In this connection it should be mentioned that an excellent review 
of complex networks in a wide range of systems in wild nature and human society was 
made in Ref. 2^. A method for embedding scale-free networks in Euclidean lattices was 
proposed in Ref. j^]. The critical dynamics, Glauber type and Kawasaki type, on two 



typical small-world networks, adding type and rewriting type, was studied in Ref. 
Newman considered assortative mixing of various types using empirical network data, 
analytic models, and numerical simulation. He showed that assortative (or disassortative) 
mixing is indeed present in many networks, it can be measured, and its effect on network 
structure and behavior can be examined. Generalized the Bethe-Peierls approach to random 
networks with degree correlations, and the analysis of the vertex cover (VC) problem as 
a prototype optimization problems defined over graphs |23l] has shown that uncorrelated 
power-law networks are simple from the point of view of combinatorial optimization, and 
inhomogeneities of neighboring vertices can be exploited. 

In hierarchical networks, the degree of clustering characterizing different groups follows 
a strict scaling law, which can be used to identify the presence of hierarchical organiza- 
tion in real networks [s^. The consideration of structured scale free- networks restores the 
order-disorder transitions in spite of the hubs, but the value of the order parameter for the 

"3^1 . A nonstandard, susceptible- 



disordered state reveals the existence of ordered clusters 
exposed-infected (SEI), compartmental model for disease epidemics comprising latency and 

;emporal decay in the rates of infection, also known as quenching, has been examined in Ref. 

3^ . Eventually, one can recognize that the study and solutions of epidemic models has been 
of a long-standing interest in theoretical physicists motivated by analogy between spreading 
phenomena in epidemic and physical systems and by application of analytical techniques. 

Numerous methods are successfully used in statistical physics to describe the distinc- 
tive characteristics of chaotic dynamics of various networks. However, three vexing features 
that are difficult for the detailed analysis can be observed in real networks. Among them: 
nonstationarity, nonlinearity, and nonequilibrium phenomena. Furthermore, the significant 
peculiarities of networks are directly related to the discreteness in time of object-subject 
registration response. Non-Markov and long-range statistical memory effects also play the 
leading part in epidemic dynamics in social networks. For this reason we offer a new sta- 
tistical approach in the research of epidemic processes with the help of the our statistical 
theory of non-Markov stochastic discrete processes in this work. This method 

has already been used by the authors to study various real problems in cardiology js^ [s^ , 
seismology [s^, neuropsychology js^ and neurophysiology [s^. With the help of this 
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method we intend to find the full spectrum of the dynamic kinetic and spectral character- 
istics of the studied complex system. We shall receive a more detailed representation of 
the stochastic dynamics of epidemic processes taking into account the discrete behavior and 
nonstationarity of the dynamic and kinetic epidemic parameters. 

This work deals with the study of statistical properties of epidemic dynamics. The key 
distinctions of grippe and respiratory infections epidemic processes are defined with the help 
of the long-range memory effects and statistical effects of non-Markovity in these processes. 
In particular, one of the key problems consists in revealing the role of randomness, regularity 
and predictability of epidemic processes. 

The paper has the following structure. In Section 2 we show the basic points of our 
statistical theory of nonstationary discrete non-Markov processes in complex systems jssi 
34 , . The description of the experimental data and the technique of their processing are 
given in Section 3. The received results and their analysis are submitted in Section 4. The 
conclusions are submitted in Section 5. 



II. STATISTICAL THEORY OF NONSTATIONARY DISCRETE NON-MARKOV 
PROCESSES IN COMPLEX SYSTEMS. BASIC CONCEPTS AND DEFINITION 

A brief description of the theory is presented here. We will consider a spread of epidemic 
as a time evolution of a patient's number (discrete random variable x{T)). Therefore we 
can use results of the statistical theory of discrete non-Markov non- stationary processes 
for complex systems. The developed description of the theory and the representation of 
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quantities used here can be found in the works published earlier 

While analyzing complex systems we obtain discrete equidistant series of experimental 
data, the so-called random variable 

X = {x{T), x{T + r), x(T + 2r), ■ ■ • , x(T + fcr), ■ ■ ■ , x{T + tN - r)}. (1) 

It corresponds to a time series of measured signal during the time (A^ — 1)^, where r is 
discretization time of the signal. In this work we take the number of patients per day x{T) 
as a measurable parameter, {N — l)r is a time interval of data recording (6 years), r is one 
day. 

For the dynamical analysis, it is more convenient to use a normalized time correlation 
function (TCF). For nonstationary discrete processes TCF has the following form {t = 
mr, N — 1 > m > 1) 

^ N-l-m 

"'"- (iv-n^MOMi) g &(r + ,.)fa(r + + (2) 
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where a{0) and a{t) is the variances of the initial (t = 0) and final (at moment t) dynamic 
states of the systems, correspondingly. The properties of TCF a{t) are determined by the 
conditions of normalization (at t=0) and attenuation of correlations (at t oo) 

lima(t) = 1, lim a(t) = 0. (3) 

t^O t-^oo 

If we take into account the non-stationarity and discreteness of complex systems for real 

brm of a closed set of the finite- 



processes, the kinetic equation for TCF a{t) has the 
difference kinetic equations of the non-Markov type 
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= Aia(t) - rAi Mr{jr)a{t - jr). (4) 
i=o 

Here Ai is a relaxation parameter with the dimension of square of frequency, and parameter 
Ai describes the eigen-spectrum of Liouville's quasioperator L 

_ . <A°(0)LA°(0)> _ < A°(0)Li2L2iA°(0) > 
<|A°(0)P> ' <|AO(0)P> ' 

where angular brackets mean a scalar product of state vectors. 

The function Mi (jr) in the rhs of Eq. represents the first order memory function 

< A°(0)L,2L2iA»(0) > • ' ' 

In Eqs. (0) and © operator L is a finite-difference operator 

A 

iL = -— , At = r. 

At' 

where r is a discretization time step, Lij = IliLIlj { i,j = 1,2) are matrix elements of 
splittable Liouville's quasioperator, Hi =11, n2 = P=l— 11 and 11 are projection operators 
(for more details, see, Refs. js^-js^l). It is easy to see, that in Eq. ^ we deal with the 
time correlation of new orthogonal dynamic variable Z/2iA°(0). 

Eq. represents the first equation in the chain of finite-difference kinetic equations 
with memory for the discrete TCF a{t). One can recognize that the memory function Mi{t) 
takes the statistical memory about previous states of the system into account. By using 
Gram-Schmidt orthogonalization procedure [3^ we receive the recurrent formula, in which 
the senior dynamic variable W„ = W„(t) is connected to the junior one in the following 
way: 

Wo = AO(0), Wi = {zL-Ai}Wo,... W„ = {^L-A„_i}W„_i+A„_iW„_2 + ..., n > 1. 

(7) 
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Introducing the corresponding projection operators we come to the following chain of con- 
nected non-Markov finite-difference kinetic equations (t = mr, n = 1,2,...) 



AM (t) 

^ = K+iMn{t) - rA„+i Mn+^{jT)Mn{t - jr). (8) 

3=0 

Here A^+i is an eigen-value of the Liouville's quasioperator, and relaxation parameters A„_|_i 
are determined as follows 

_ . < W„LW„> < W„_i(2L-A„)W„> 

<|w„|2> ' <|W„_i|2> 

In this work we analyze the short time series with the help of Eqs. (^-(jH)) where non- 
stationary functions may not be taken into account js^. So, we shall use full set of the 
dynamic, kinetic and spectral parameters, functions and characteristics as a quantitative 
information measure to describe of spread of epidemic in human population. Among them, 
here are temporal dynamics of the first four orthogonal dynamic variables, phase portraits 
in some planes of dynamic variables, TCP, junior memory functions and their power spec- 
tra, frequency dependence of first three points of non-Markovity parameter, and local time 
behavior of the locally average relaxation parameters. 



III. EXPERIMENTAL DATA AND DATA PROCESSING 



The registration and accounting materials, as well as the results of certain studies are the 
initial data for epidemic research in medicine. We have made use of the accounting materials 
of Kazan sanitary - epidemiological stations located in industrial districts. The concrete data 
were given by the sanitary - epidemiological station of the Privolzhskii district of Kazan 
covering the period from 10.27.1995 to 03.05.2002. Figs. l.A presents the experimental raw 
data of epidemic spread, involving the total population of the Privolzhskii district (about 
300 thousand people). The first group of the data represents a six-year dynamics of grippe, 
the second group describes a six-year dynamics of acute respiratory track infections in this 
district. The obtained data were processed with the help of the above introduced technique. 
The set of three memory functions was calculated for each sequence of the data. The power 
spectra for each of these functions are obtained by the fast Fourier transform (FFT). We 
will also show the phase portraits in plane projections of the multidimensional space of the 
dynamic orthogonal variables. For a more detailed diagnosis of the system we will consider 
the frequency spectrum of the first three points of the statistical spectrum of the non- 
Markovity parameter. In this study we will use the frequency dependence of the statistical 
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spectrum of non-Markovity parameter 
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(9) 



Here i = 1,2.. and fii{uj) is a power spectrum for ith order memory function of the ith 
relaxation level. The statistical spectrum of non-Markovity parameter ei{uj) constitutes an 
information measure of Markovity and non-Markovity or randomness and regularity in time 
evolution of underlying systems. 

IV. DISCUSSION OF THE RESULTS 

In this section the quantitative and comparative analysis of the six-year period of grippe 
chaotic dynamics (Fig. l.B.a) and of ARTI (Fig. l.B.e) will be carried out on the basis of the 
theory submitted in Section 3. The time series of the initial signals, the phase portraits of 
the dynamic variables, the power spectra of TCF and the junior functions of memory as well 
as frequency dependence of the first three points of the statistical non-Markov parameter 
are submitted in Figs. 1-4. We develop a new approach in the study of epidemic processes 
on the basis the local time behavior of the relaxation and kinetic parameters Ai, A2, A3, Ai 
and A2 (Figs. 5, 6) also. 

The six-year sampling of grippe and ARTI dynamics are submitted in Fig. l.B. (Figs. I.E. 
a-d for grippe. Figs. l.B. e-h for ARTI). The time series of the orthogonal variable Wq (Fig. 
l.B.a), Wi (Fig. l.B.b), W2 (Fig. l.B.c), W3, (Fig. l.B.d) for the chaotic dynamics of grippe 
and ARTI have appreciable symmetry relative to straight line Wi = 0, where i = 1...3. The 
analysis of the time series of variable Wo (Fig. l.B.a) for a grippe epidemic shows that the 
moment of the grippe outburst or appearance of epidemic corresponds to the most significant 
fluctuations of this dynamic variable. Epidemic outbursts of the disease are located at equal 
distances from each other. This is the evidence of periodicity and interval constancy between 
grippe epidemics within the whole time of observation. The time series for three orthogonal 
variables Wi (Fig. l.B.b), W2 (Fig. l.B.c), W3 (Fig. l.B.d) of grippe epidemics are almost 
perfectly symmetrical. The time series of dynamic variables PFq (Fig. l.B.e), W^i (Fig. l.B.f), 
W2 (Fig. l.B.g), W3, (Fig. l.B.h) of ARTI epidemiological processes present a different 
picture. The visible fluctuations of the dynamic variables for ARTI are observed in the 
autumn-winter period same as in case of grippe epidemics. However, dynamic noises can be 
observed between the outbursts of ARTI epidemics. They are represented by monotonously 
inequable parabolas. One of the half-parabola corresponds to one of the epidemics and 
another corresponds to the subsequent one. The time development of orthogonal variables 
Wi (Fig. l.B.f), W2 (Fig. l.B.g), W3 (Fig. l.B.h) of the ARTI epidemiological process 
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represents symmetric series about straight line Wi — 0. There are seven outbursts in the time 
behaviour of the dynamic variables in the seven autumn - winter periods (from 10.27.1995 
to 03.05.02). They corresponds to the seven cycles of grippe and ARTI epidemics. However 
the time series of orthogonal variables Wi {i — 1,2,3) for ARTI differ from grippe's case by 
the strongly expressed asymmetry and the existence of a clearly expressed noise. 

In Fig. 2 the phase clouds in the six plane projections of the four first dynamic variables 
Wi, i — 0...3 for epidemic processes of grippe and ARTI are submitted. In case of grippe the 
phase clouds have a well defined asymmetry about the centre of coordinates. All the phase 
clouds contain a centralized nucleus and an aggregate of points scattered on perimeter in a 
fan-shaped way. The interval of the dispersal makes up 900 r. Owing to such asymmetric 
stratification the phase clouds get an elongated shape. The phase portraits of the ARTI 
epidemiological process have smaller symmetry about the centre of coordinates. The nucleus 
of the portrait gets the shape of a right angle and resembles a " boot" . Apical emissions of 
separate points are appreciable on each side of the "boot". The interval of the dispersal 
of the points decreases up to 700 r. The stratification of the phase clouds of a grippe 
epidemiological process is more essential, than in case of an ARTI epidemic. The analysis of 
the data in other Kazan districts reveals a prominent phase portrait of a grippe epidemic. It 
is the elongation of the phase portrait for an ARTI epidemic namely its "boot- type" shape. 

In Fig. 3 the power spectra of TCF /Xo(<^) and three junior memory functions ^i{uj),i — 
1, 2, 3 for grippe and ARTI epidemiological processes are submitted. The frequency spectra 
of /io(<^) are given in doubly-log scale for a more detailed analysis of the data. In case of 
grippe the fractal dependence in the region of intermediate and high frequencies is found in 
the power spectrum of the initial TCF ^^{oj). The area of these frequencies is divided by a 
small outburst of power on frequency a; = 3 * 10~^/.w. {If.u. — 2t:/t). In the domain of the 
low frequencies sharp breaks of the spectrum are observed. The frequency dependency of the 
next three junior memory functions iii{u!),i — 1, 2, 3 are accompanied by spectral outbursts 
at equal intervals and by condensation of spectral lines near the outbursts. A spectral noise 
have been presented in all diagrams of the frequency dependencies. The spectrum of the 
initial TCF //o(<^) for an ARTI epidemic process differs from a grippe epidemic only in the 
high frequency area. Additional spectral peaks are appreciable at frequencies cu — 10~^f.u. 
and uj — 10~^^^f.u. The power spectra of the three junior memory functions fj,i{uj),i = 1,2,3 
for an ARTI epidemic have bright fractal dependence with several spectral peaks at every 
0.15/.-^. 

In Fig. 4 the spectra of the first three points of statistical non-Markov parameter 
where i — 1,2,3 are presented. Fractalness is appreciable only in the behaviour of 
With approximation to frequency a; = the monotonous amplification of markovization is 
observed. The spectra of the non-Markovity parameter for the second point of 62(00) for 
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ARTI and grippe epidemics accept value £2(<^) = 1 are almost identical but it have small 
fluctuations. It testifies to strong non-Markovity. The sinuosity with two frequency crests 
in area a; = 0.1 f.u. and uj — 0.3 f.u. is appreciable in frequency spectra of the third point 
of statistical non-Markovity parameter £3(0;). 

In Fig. 5 the time dependence of local (with local averaging) kinetic and relaxation 
parameters Ai, A2, A3, Ai and A2 for a grippe epidemic is submitted. The sampling for 
the local averaging includes 200 points. The similar size of the samphng allows to consider 
the physical nature of a real object (a grippe and ARTI epidemics) in a more detailed way. 
Breaks are observed at a great number of points in the time scanning of the parameters. 
In the time dependence of all the parameters within six years of every autumn - winter 
period the system drops out of the steady state and returns to it. The whole process can be 
presented as recurrence of two various processes. At first small fluctuations of the parameter 
are observed. The fluctuations visibly grow in the moment of seasonal epidemics. Then the 
system acquires the state of quasistable balance. The time interval for each process is 125 
points, that corresponds precisely to one of the autumn - winter periods. 

Recently the correlation analysis has experienced a marked lack of information concerning 
the object of the research. Procedure of local averaging of various parameters allows to 
examine the separate hidden properties of studied objects. The characteristic feature of the 
usual correlation analysis is the fact that the greatest possible set of signals is required for 
the qualitative analysis of the properties of the object of the research. At longer sample of 
such signals it is possible to receive more exact information with the help of the correlation 
analysis. Let us take a random non-Markov process as an example. This process consists 
of a sequence of alternating states. Thus the problem of extraction more information not 
only about the common process but about various separate states of the system arises. In 
this case the use of the correlation analysis for all the time series will be inefficient. The 
processing of the signals is necessary for separate local sites of the whole time series. It will 
allow to consider the properties of separate dynamic states of the system. 

Hereinafter a new method of data processing based on the local averaging of kinetic and 
relaxation parameters is offered. This method allows to consider the properties of separate, 
non-stationary states of the systems. The idea of a method is the following: there exists 
an initial data set. Let's take a samphng in length N of signals and calculate its kinetic 
and relaxation parameters. Then the operation of " step- by-step shift to the right" for one 
time interval is carried out. The kinetic and relaxation parameters are calculated again. 
The "step-by-step shift to the right" is continued to the end of a time series. Such locally 
averaged parameters have high sensitivity to the effects of intermittency and nonstationarity. 
If the initial time series has some irregularity, it is instantly reflected in the behavior of the 
locally averaged parameters. 
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The use of this method requires the choice of the optimal length of a sampling which 
enables to receive the most trustworthy information. If a sampling is too short, noise effects 
do not allow to receive qualitative information. Besides with a short length sampling we 
have significant errors. On the other hand at a great length of a sampling locally averaged 
parameters lose " sensitivity" necessary for the study. As a result of the studies of different 
lengths of local sampling we have received the optimal length which makes 100-120 points. 
Further proofs of all aforesaid will be given below. 

In Fig. 6 the time dependence of local kinetic and relaxation parameters Ai, A2, A3, Ai 
and A2 for an ARTI epidemic is submitted. 200 points where selected for the procedure 
of localization. Due to a more significant noise of the ARTI chaotic dynamics the time 
scanning of then parameters essentially differs from grippe. The time dependence of the first 
parameter Ai refiects the physical features of an ARTI epidemic process. Instead of small 
fiuctuations at the beginning of the epidemic sharp fiuctuation splashes appear. Between 
them the system acquires a rather quiet state followed by a new epidemic peak. When the 
epidemic is over only small fiuctuations in the system can be appreciable. This process 
repeats periodically during all long-term interval. Some definite regularity can be observed 
in the time dependence of the next kinetic and relaxation parameters. The sharp peak at 
the beginning and at the end of the ARTI epidemic process relates to them. Only small 
fiuctuations of the system are appreciable at the time of epidemic. After each ARTI epidemic 
the considered system comes back to the state of relative balance. 

Judging by the behavior of the parameters the distinction between the epidemic of grippe 
(I) and that of ARTI (II) comes to the following. The initial values of the intensity //^(a; = 0) 
for alH = 0, 1, 2, 3 for systems I and II differ by more than one order. It is the evidence 
of remarkable distinction in times of relaxation (correlations). For system I these times are 
almost 20-50 times shorter, than for system II. The intensity scales also differ almost by 50 
times. Thus the initial values of £1(0; = 0) differ almost by 5 times for systems I and II. 

On the basis of the aforesaid it is possible to come to certain conclusions about a great 
randomness of grippe dynamics in comparison with ARTI dynamics. The last one has greater 
regularity and predictability. Other feature of the investigated epidemics is the existence of 
the three additional well - defined spectral splashes except the basic annual. From figs. 3 
b-d and 3 f-h one can see that these groups of spectral peaks correspond to the characteristic 
frequencies equal accordingly to 0.14 f.u., 0.2Sf.u. and 0.43 f.u. The intensity of the high- 
frequency peak is almost an order lower, than the two other peaks of lower frequencies. 

The strong distinction in the dynamic behavior of grippe and ARTI epidemics is shown 
in local time dependence of relaxation constants Ai , A2, A3, Ai and A2. It is necessary to 
remind, that parameters Aj, i = 1,2,3... remind Luapunov's exponents for various relaxation 
levels. For system II great regularity is observed. All parameters of Ai, A2, A3 are negative. 
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which is the evidence of the relative stabihty of the system. For system I (grippe) owing to 
the general instability the parameter Ai(t) gets small positive values on tail. It corresponds 
to the slump of epidemic. From Eq. (8) one can see, the system becomes unstable at values 
Ai(t) > 0. 

The time behavior of parameters Ai(t) and A2(t) is more dramatic. If Aj(t) > 0,2 = 1,2 
then solution of Eq. (8) corresponds to a steady condition of social network, if Aj(t) < 
instability grows. From Figs. 5 d, e and 6 d, e it is clear, that annual epidemics of grippe 
have a greater instability of the system. In case of ARTI the behavior of A2(t) everywhere 
has steady state of system II. However the parameter of Ai(t) in the greater part of the time 
scale is negative, this testifies to a loss of stability in this area. 

So, according to our theory, parameter Ai(t)= R make sense a relaxation rate of random 
process. In case of grippe (I) isoline exists in behavior Ai=Ai(t) with values 0.08 r^^ < -R < 
0.24r~^. Simultaneously, a relaxation rate on recession of grippe epidemic sharply grows up 
to (0.8-1.0)r~^, that it is in 4 - 12 times, approximately. This implies the greater stabilization 
of system, than during rise or on a peak of epidemic (see, fig. 7 a, b for more details). 

In case ARTI isoline in behavior Ai(t) does not exist at all. It testifies about quasist ability 
and seasonal prevalence in relaxation behavior, in - whole. Absolute values of relaxation 
rates fluctuate in an interval 0.24 r^^ < -R < 0.6r^^. On a tail a relaxation rate on recession 
of epidemic ARTI speed of grows in 2-2.5 times. Therefore system, as a whole, has the 
greater macroscopic stability. 

We can draw the following conclusion from comparison of relaxation parameter Xi{t) of 
grippe and ARTI (see, table 1 for the whole sampling). Relaxation rates differ almost in 3.4 
times. It testifies to macroscopic stability of epidemiological process of ARTI. On the other 
hand, local amplification of relaxation rate (on a tail of epidemic) in case of a grippe creates 
an additional source of stabilization for unstable, in the whole, process of grippe epidemic. 



V. CONCLUSIONS 



In this paper we have shown that the chaotic dynamics of epidemic processes of grippe and 
ARTI in real social network can be considered as a stochastic discrete non-Markov process. 
It is more convenient to use the statistical theory of the discrete non-Markov processes for 
real objects and alive systems js^ - [s^ to study the similar processes. The used theory 
allows to define essential distinctions between the epidemic processes of grippe and acute 
respiratory track infections by degrees of randomness, regularity and predictability. By 
processing the experimental data about patient's numbers we received the time series of 
dynamic variables Wi{t) and calculated the memory functions and the first three points of 
the statistical spectra of non-Markovity parameter. We used the power spectra, which were 
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received with help of the fast Fourier transform (FFT) to analyze the diverse time functions 
(correlation and memory functions). 

The new qualitative approach in studying statistical properties of epidemic processes of 
grippe and ARTI originates while using locally average kinetic and relaxation parameters 
Ai, A2, A3, Ai and A2. 

It is possible to come to the following conclusions on the basis of the conceptions about 
long-range memory and localization of parameters. The epidemic process of grippe causes 
a great danger to mankind because of its randomness and unpredictability. Each epidemic 
process of grippe is characterized by a sudden outburst and sudden attenuation. The peak 
of an epidemic lasts a comparatively short period of time (60 days). It means efforts to 
liquidate the consequences of a grippe epidemic should be made during this period of time. 
The epidemic process of ARTI is characterized by a greater regularity and predictabihty. 
It is accompanied by a small outburst at the beginning of the epidemic, then the number 
of patients diminishes. After it a drastic outburst corresponding to the epidemic peak is 
observed. The process ends with small short-term splashes which later quickly fade. Unlike 
grippe, an ARTI epidemic is characterized by a less drastic outburst and attenuation. During 
the underlying interval insignificant outbursts of respiratory infections are observed followed 
by periodic outbursts of ARTI. 

We have revealed and explained essential distinctions in phase portraits, frequency spec- 
tra and quantitative estimations of memory functions and statistical non-Markov parameters 
between the epidemic processes of grippe and ARTI. Relative distinctions of the quantitative 
data about epidemics of grippe and ARTI are the evidence of the distinction in long-range 
memory, order and organization of each process in the certain social system. These distinc- 
tions allow to define the degree of importance of each of the epidemics. In particular, we 
have shown that the epidemic process of grippe and respiratory infections have essentially 
different quantitative measure of randomness and regularity. 

The local kinetic and relaxation parameters of epidemic processes of grippe and ARTI 
allow to study statistical features of these systems in detail. The local time dependencies 
allow to find the internal features of the epidemic process. It helps in the study of real 
objects. With abundant experimental data the developed method allows to define the laws 
of epidemic processes. 

The received results can be of practical value while studying other diseases and compar- 
ative estimation of epidemic danger. 

In this paper we have demonstrated that the set of relaxation, kinetic and spectral pa- 
rameters and characteristics of discrete non-Markov stochastic processes are valuable for the 
description of the role of randomness, regularity and predictability of epidemic processes, 
the spread of grippe and ARTI in the world. 
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Since the similar situation is typical of the majority of epidemic diseases on networks 
our conclusions are of profound importance for a large number of physical, biological and 
technological networks. The results received here are connected with non-stationary and non- 
ergodic processes in stochastic systems. Therefore they are certainly of value for physics of 
the disorder matter in which similar processes can occur on molecular or structural levels. 

Distribution of atypical pneumonia (SARS) attaches special value to any researches of 
ARTI. The cases similar SARS, unfortunately can a multiply in the future . In connection 
with stated above a study of ARTI infections represents the important scientific significance. 
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VII. FIGURES CAPTION 

FIG. l.A. Six-year temporal dynamics of the experimental raw data of grippe (l.A.a) 
and acute respiratory track infections (l.A.b) in Privolzhskii district of Kazan, Russia, from 
10.27.95 to 03.05.02. The dynamics of ARTI is initially characterized by a slow increase and 
then by slow decrease in the number of patients. With grippe we have a more sharp picture, 
characterized by low level of the background noise, then by a sudden increase and decrease 
of the disease. 

FIG. l.B. Time series of orthogonal variables Wq (a, e), Wi (b, f), W2 (c, g), W3 (d, h) of 
chaotic dynamics of grippe (a, b, c and d) and acute respiratory track infections (ARTI)(e,f,g 
and h) . The spectral peaks in time series display the seasonal periodic occurrence of grippe 
and ARTI epidemics. The time series of orthogonal variables for the epidemic process of 
ARTI differ due to the occurrence of an appreciable noise. 

FIG. 2. The phase clouds for epidemic processes of grippe and acute respiratory track 
infections (ARTI) in plane projections for the various combinations of orthogonal variables 
Wi, Wj where i,j — 0, 1...3. During the epidemic of grippe the phase points disperse about 
the nucleus. The phase clouds for the epidemic process of ARTI acquires the shape of a 
"boot" and are grouped about nucleus. 

FIG. 3. The power spectra of time correlation function (TCF) /xo(a;)(a) and of the three 
junior memory functions ^i{uj),i — 1(6), 2(c), 3 (d) for the epidemic processes of grippe and 
acute respiratory track infections (ARTI) . For a more detailed analysis the power spectrum 
of TCF /xo(a;)(a) is submitted in double logarithmic scale. In the field of high frequencies 
three additional spectral splashes are observed. The spectral peaks in the registered spectra 
are the evidence of of long-range memory in the system. 

FIG. 4. The frequency dependence of the first three points of statistical non-Markov pa- 
rameter £i{uj) for epidemic processes of grippe and acute respiratory track infections (ARTI). 
For the first point of £1(0;) the splash of markovization process is observe in the vicinity of 
value uj — The values of other parameters e2{oj) and £3(0;) change in a narrow interval 
about a unit. This reflects the existence of long-range memory order in the system. 

FIG. 5. The time dependence of local kinetic and relaxation parameters Ai, A2, A3, Ai 
and A2 for the epidemic process of grippe in units of r. The sampling of local averaging 
is equal to 200 points. The sudden dropout of the system from a balanced condition is 
determined by the epidemic of grippe. The reversion of the system into a former steady 
state defines the intervals between the epidemics of grippe. 
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FIG. 6. The time dependence of local kinetic and relaxation parameters Ai, A2, A3, Ai 
and A2 for epidemic process of acute respiratory track infections (ARTI) in units of r. The 
length of the sampling of local averaging is equal to 200 points. The sharp splash in a time 
series demonstrates the beginning of an ARTI epidemic. The end of the epidemic process 
is defined by the following peak. At the moment of epidemic the system is in an unstable 
state. 

FIG. 7. Time dependence of local relaxation rate Xi{t) for grippe (b) and ARTI (d) with 
initial row data for grippe (a) and ARTI (c). We show extracted sampling for comparison. 
One can see the dramatic revising of Xi{t) during both epidemic. Time behavior Xi{t) at 
grippe contains isoline's regime whereas one can notice continuous alteration of Ai (t) in case 
of ARTI. 
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